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Dissipative effects can lead to a friction term in tfie equation of motion for an infiaton field during 
the inflationary era. The friction term may be linear and localised, in which case it is described by 
a dissipation coefficient. The dissipation coefficient is calculated here in a supersymmetric model 
with a two stage decay process in which the infiaton decays into a thermal gas of light particles 
through a heavy intermediate. At low temperatures, the dissipation coefficient oc in a thermal 
approximation. Results are also given for a non-equilibrium anzatz. The dissipation coefficient is 
consistent with a warm inflationary regime for moderate (~ O.f) values of the coupling constants. 

o 
o 

(N 

, ^ , Inflationary models have made a huge contribution to our understanding of the early universe and the origin of the 
basic cosmological features which are observed today in ever increasing detail P, 0, U . This paper is concerned with 
the mechanism for converting the vacuum energy which dominates the universe during inflation into the radiation 
which dominates after inflation. Radiation can be produced during inflation, at the end of inflation or after inflation 
has finished. We will consider the first possibihty in a particular class of supersymmetric models, to see whether a 
I significant amount of radiation can be production during infiation. When this happens, the situation is described as 
*>. warm inflation HQ (see also @). 

I The production of radiation is balanced by dissipative terms in the equation of motion of the homogeneous inflaton 
^ ■ field (j). The simplest possibility is a linear dissipation term, 

(N ■ 

g: 0+(3i7 + r)0 + K^ = (1) 

] where is the (f> derivative of the inflaton potential, H is the expansion rate of the universe and F, which in general 

■ will depend on both cj) and temperature T, is called the dissipation coefficient. The dissipation coefficient was originally 
introduced as an arbitrary parameter into discussions of reheating '7| . Explicit formulae for the dissipation coefficient 
were obtained from analysing particle production j^,^] and linear response theory .10.1 . Later work has tended to use 

O • the Schwinger-Keldeysh approach which we use here, eg jUIll, El, Q EE EE HE US ■ 

D , Dissipation can lead to important consequences when T > H (called the weak regime of warm infiation) and 

■ when T > H (the strong regime) . In the weak regime, the primordial density perturbation spectrum is determined by 
^ [ thermal fluctuations rather than vacuum fluctuations. In the strong regime, restrictions on the gradient of the inflaton 

• '~j ' potential may be relaxed Q . This is especially important in the context of i^-term supersymmetric inflationary models 

rS ' (see [Tgj i for a review). The potentials in these models contain supergravity corrections which rule out the ordinary 

^ . type of inflation, but not warm inflation [T^ [20j . 

■ - - ' On the other hand, the calculations which have been done so far have assumed that radiation remains close to 
thermal equilibrium, and this leads to additional consistency requirements: 

• Thermal corrections to the inflaton potential have to be consistent with the restrictions on the gradient of the 
potential. 

• The thermalisation time r of the radiation should be smaller than the expansion timescale of the universe. 

These conditions lead to strong restrictions on the existence of warm inflation. Nevertheless, it is useful to investigate 
the close to equilibrium regime before moving on to the more difficult far from equilibrium regime. 

The transfer of energy from the infiaton to another field x is dependent the coupling strength of the interaction, 
the relative sizes of the mass of the inflaton m^, the mass of decay product and the temperature T. The theory 
of reheating is well-developed for the post-infiationary era when the inflaton is oscillating and the parameter range 
7710 > 771^. According to linear theory, the rate of decay of these oscillations is determined by the decay rate for 
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processes such as (/> ^ 2% 0, 0, 11^ • The theory of preheatine is the modification which takes into account 
non-hnear effects close to the particle production thresholds l25l l2a | . 

The theory of the dissipation coefhcient in the inflationary era is also quite well-developed, especially at high 
temperature, when T > m^. This is the case in which the consistency requirements for an inflationary model are 
most severe, with large thermal corrections when the radiation is strongly coupled and large relaxation times when 
the radiation is weakly coupled ^27] . 

Recently, there has been interest in the possibility that the inflaton can dissipate energy at low temperature via a 
two-stage decay 4> —> X ~* J/y? with the heavy field x acting as an intermediary in producing a light field y 28, 29, 30j. 
The dissipation in these models can become a significant feature if the coupling strengths of the inflaton to the heavy 
field is moderately large. 

Inflation is usually associated with a very weakly coupled inflaton to protect the gradient of the potential from 
picking up large quantum corrections. However, in supersymmetric models the restrictions on the couplings are less 
severe. The gradient of the inflaton potential can be small enough to allow inflation whilst the thermalisation time 
for the radiation remains short, allowing warm inflation for moderate values of the coupling constants in these models 

[23,E1IS2. 

In this paper, we shall calculate the dissipation coefficients in supersymmetric models which have an inflaton 
together with multiplets of heavy and light fields. Different contributions to the final result can be applied also to 
non-supersymmetric models. 

In the high t emp erature regime T > m^, our results are similar in many ways to those obtained previously for 
self- interacting [13 and coupled systems [33 ■ There are some differences for T ~ due to the details of the x^ VV 
decay. These differences affect the temperature dependence of the dissipation coefficient. This could potentially have 
observable consequences, since studies of density perturbations in models with different forms of dissipation terms have 
shown that the spectrum of density perturbations is quite sensitive to the temperature dependence of the dissipation 
coefficient 0. 

We have new results for the regime H < T < m-^ and < m^. Previous calculations in this regime have been, 
in effect, highly non-thermal [ME9- Our results differ substantially from these non-thermal calculations. We find 
that the leading order dissipation terms in the infiaton equation vanish in the low temperature limit, but they are 
sufficiently large to ensure that the temperature does not fall to zero during inflation. We shall see in section VI that 
warm inflation appears to be allowed for moderate values of the coupling constants. 



II. THERMAL FIELD THEORY 



We begin with some of the basic notions and results from thermal field theory which will be used later on in this 
paper. More details and proofs can be found in standard texts, eg |35| . 
The ensemble average of a Heisenberg operator A can be written 

^ = ^s:<|}|!M = '-^ (2) 

Pii'iVim'} tr(p) 

where p is the density matrix and repeated indices indicate sums. We shall be using a shifted version of the Schwinger- 
Keldysh, or closed-time path formalism to evaluate ensemble averages [3^[33- For a scalar field cj), the generating 
functional 



Z[(j)i,(j)2,Ji,J2]=Pii'{ipi\T*eyip(^-i J J2{(f) ~ (j)2)^ T exp (^i J Ji{(j) - 



1^.'), (3) 



where T denotes time ordering of the operators with the smallest time on the right and T* denotes time ordering of 
the operators with the smallest time on the left. Shifting the field operators by (pi and (j)2 is done for later convenience. 

Following Calzetta and Hu fs^l, we combine the sources into a pair Ja = (Ji, J2) and introduce a metric Cab = 
diag(l, —1) to raise and lower the index a. This removes the inconvenience with the relative signs of the two source 
terms. The generating functional gives rise to four two-point functions. 



^ab — ^ 



Z SJ'^SJ'' 



(4) 

(Jl,J2) = 



If we let f] — (j) — (j)a, these two-point functions can be expressed in terms of the functions 

G>{x,x') = G<{x',x) = i{fi{x)fi{x')). (5) 
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The individual components are 

Gii{x,x') = e{t~t')G>{x,x') +e{t' -t)G<{x,x') (6) 

G22{x,x') = e{t~t')G<{x,x') + 0it' ~t)G>{x,x') (7) 

Gi2{x,x') = G<{x,x') (8) 

G2i{x,x') - G>ix,x') (9) 

The two-point function Gn is the thermal analogue of the Feynman green function Gp- (Placing a factor of i in front 
of the expectation values in the definition of the two-point functions follows 39] but differs from most references.) 

Consider a system which remains close to local equilibrium. The natural starting point for perturbation theory is 
a free theory with an initially thermal distribution of states. The two-point functions for the interacting theory are 
given by Feynman diagram expansions with lines representing free theory propagators —iGoab and vertices derived 
from an interaction Lagrangian. The mass of the propagators and the vertices depend on the background fields 0a, 
which may be time or space dependent. Each of these vertices has a label a attached, and the final expression for the 
diagram has a sum over the a's. Later, we will make use of the following result, 

• The maximum time rule. The contribution from a diagram vanishes if the time at any vertex is larger than all 
of the times on the external lines of the diagram. 

This follows directly from the operator form of the two-point functions H6I9|1 . 

The next step is to consider the background fields (j)a- These can be identified with the ensemble averages of the 
field operator by setting 

We assume that this equation can be solved for Ja[4>b]- The effective field equation is then F = 0, where 



Fmx) 



0, (11) 



S4ibix) 
and the effective action 

^[^a] = -^lrlZ[cf>a,Jbm]■ (12) 

The effective action can be expanded, as usual, by a sum over 1-particle irreducible vacuum diagrams, with internal 
lines representing the free-theory propagators —iGoab- 

Dirac and Majorana fields can be included in the Schwinger-Kcldysh formalism without any difficulty. We define 
the fermion two-point functions 

S>{x,x') = i(V'(x)V^(x')), S<{x,x') = -i(^{x')i}{x)) (13) 

and 

Sii{x,x') = 9{t-t')S>{x,x') +0{t' -t)S<{x,x') (14) 

S22ix,x') = 9{t ~ t') S<{x,x') + 9{t' - t) S>{x,x') (15) 

Si2ix,x') = S<{x,x') (16) 

S2iix,x') = S>{x,x') (17) 

where S^{x,x') = -S>{x',x) = S<{x,x'). (For a Dirac matrix Af, let M = 7°Af^7° and M = {C-^MC)^, where G 
is the charge conjugation matrix; ^ = —7^.) 

Finally, we consider an interacting system in thermal equilibrium with temperature T. The background fields must 
then be constant, = (f>2 and the scalar two-point functions Gab{x,x') = Gab{x — x'). We set 

Gab{P) = j d^xGab{x,x')e-'P<---'^ (18) 

where the 4-momentum P — (p,a;). The vacuum polarization or self-energy Hab is defined by 

G-l°''(P) = (P2 +TO2)^afc ^ jjab^ (19) 

A diagrammatic representation for the two-point functions is shown in figure ^ 
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FIG. 1: Diagrammatic representation for the free two-point function —iGoab, the two-point function —iGab and the vacuum 
energy -illab- 



We shall often represent the dissipative part of the self-energy by a function a, defined by 

a-i(n2i-ni2) (20) 

The function a is related to the relaxation time Tp of the particle state via the relation a — 2ujpTp^, where ujp is the 
energy of a state with momentum p. 

In the case of fermions, we denote the self-energy by Safe, and 

S-^^\p) ^ . p + m)c'''' + i:''^ (21) 

Further properties of the equilibrium propagators are listed in the appendix. 

III. DISSIPATION TERMS 

A. Adiabatic approximations 

We shall examine the effective field equation for an inflaton field which is interacting with radiation and evolving 
very slowly compared to other dynamical timescales. We shall construct an adiabatic approximation in the form of a 
derivative expansion. 

Some conditions will be imposed on the radiation. The thermalisation timescale which governs the return to 
equilibrium of a small non-thermal fluctuation in the radiation will be denoted by , the evolution timescale of the 
inflaton field will denoted by and the expansion timescale of the universe by r^. We take 

• < T0 

• < Tq 

• r> ij 

The first two conditions are necessary if the system is to remain close to thermal equilibrium. The second and third 
conditions allow us to neglect the expansion of the universe when applying the thermal field theory. Techniques exist 
for taking into account both expansion and thermal effects (see ,29, 40] for example), but typically we expect either 
the thermal or expansion effects to predominate, and in this paper we consider the former only. 
Consider the effective field equation. 



sr 



(22) 



where </> is spatially homogeneous and varies slowly about its value (j)t at time t. We can set (500 = (f'a — (f't and expand 
Fby 

CO 

F{x) = Y,Fr.{x) (23) 

n=0 

where 



1 f J""*" 

F„{x) ^ / d'^xi . . . d'^Xn 

n\ J d(/)i(x)d(/)Qi(a;i 



baAxi)---&(t>aAXn) (24) 
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The first term _Fo represents the part of the field equations which contains no derivative terms, and can be expressed 
as the derivative of an effective potential which we denote by V^. 

The next term Fi contains the derivative terms from the classical action and the equilibrium self-energy n^j,(k, t—ti) 
of the inflaton perturbations rj. According to the maximum time rule, we only need to integrate over times ti < t. 
The spatial integral is equivalent to evaluating the self-energy at zero spatial momentum, 

Fi(a;) = + 3iJ(/.- / dt, n1^iO,t - h)Sr(h)\^^^^ (25) 

J —oo 

We shall assume that the self-energy introduces a response timescale r. If (p is slowly varying on the response timescale 
r, then we can use a simple Taylor expansion and write 

= (Pit) + {h ~ t)<p{t) + ... (26) 

The inflaton equation of motion including the linear dissipative terms is then 

4>+{3H + T)<}> + V4, = (27) 

with dissipation coefficient 

/•oo 

r = 2/ dt' Re{n:^^{0, t'))t', (28) 
Jo 

where we have used = —112]^. 

This result is typical of linear response theory j41] . The range of applicability depends on the assumption that 
the system remains close to equilibrium. We also require t < for the existence of a local dissipation term. If this 
inequality breaks down, then we can replace the dissipation coefficient by the non-local expression (|25|) . although the 
non-local expression can break down when the leading term Fi is not a uniform approximation over large timescales 



B. Simple example 

We shall illustrate the adiabatic approximation with a simple example which includes a massive boson field x ^-^d 
a massless fermion ij: with an interaction Lagrangian 

£/(0, X, i') - - -^hxi^^ (29) 

The inflaton field (p is taken to be a slowly varying function of time. 



^ \/ 

V /• s. 

-2ig'^ct>a -2ig^ 
FIG. 2: Vertices for interactions between the fields (f> and x. 



The Feynman diagram vertices are generated in the shifted Schwinger-Keldeysh formalism by 

^'l{(l)a,Va,Xa,tpa) = Cl{<j)i +?7l,Xl,V'l) -^l{<l>2 +'72,X2,V'2) (30) 

The vertices for interactions between 77 and x s-re shown in figure |2| We are interested in the situation where x as a 
heavy field, with mass ~ g4>. We shall therefore regard the first vertex as 0{g). 
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2\ 



1 ; 



"Tr 



FIG. 3: Contributions to the rj self-energy of order (left) and (right). 



The dissipation coefficient can be obtained from eq. (|28|) . The contribution to the self-energy of the rj field at order 
is given by the ffist diagram in figure |31 with two x propagators, 

r = 4/02 Im y dt' iG^At')ft'. (31) 

This can also be written in terms of the spectral function p^, using (|A8D . 

/OO 7 
^^z(l-fnKe— (32) 



After applying the identity 



we get 



Im / dt't'e'"^'' = -ttS'{uj), (33) 



V 94,2/' d^k [^dw ^ , 

There are two contrasting limits in which the result can be replaced by an approximation. First there is the weakly 
coupled limit with small h and fixed T. The energy integral is dominated by the point lo — LOk, 

c.,. -(fc2+m2)V2, (35) 

which lies close to two poles in the spectral function l|A10p . These two poles are at = Wfe ± ir^^ (possibly on a 
different Riemann sheet due to branch cuts), where is the relaxation time for the y boson. The integrand can be 
expanded about = Wfc to obtain the result first derived by Hosoya and Sakagami [lOj , 

/(Pk T 
(2^2if"(" + ^)- 

The contribution to due to a single fermion loop is of order h~^, as we shall see in a later section, so that 
r = 0(.g2/i-2). 

The low temperature limit Pm^ —> 00 for fixed coupling h gives a quite different approximation. In this case, 
exponential terms control the energy and momentum integrals and the dominant contributions come from uj <C 
and k <C m^. According to equation (|A10|) . the spectral function can be replaced by 

Px « ^ (37) 

where is the low-energy decay width of the x boson. The dissipation coefficient given by equation (|31|l becomes 

r-^/?7#i32/3^F>(n-fl) (38) 



m6 J 27r (27r)3 ' ^ 

The magnitude of this term depends on the behaviour of the vacuum polarisation at small uj. We shall give results 
for in a later section. The fermion contribution to the vacuum polarisation gives a contribution F;^ ~ h^K'^m~^ , 
and F cx T^. The dominant contribution to F at low temperature comes from light scalar interactions, which we shall 
consider in section III. 
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C. Zero-temperature dissipation terms 



The simple example of the previous section suggests that the dissipation coefficient in the inflaton equation vanishes 
at order in the zero temperature limit for a non-expanding universe. This contradicts some previous results 

I28L l3Clj . but it is consistent with the conclusions of Boyanovsy et al. who found that their zero-temperature 
dissipation term could not be localised. We shall therefore examine this issue further. 

Consider the simple model again, with the 0{g^) contribution to the dissipation coefficient 

r = 'lg*(j)'lm dt'G^^it'ft' (39) 



(27r)3 

In the vacuum, G21 is the Wightman function, and includes corrections due to interactions with light fields. The 
spectral representation of the Wightman function is non-vanishing only when the energy lies in the physical particle 
spectrum, 

G2i(c^) =i0(c^)PxM- (40) 

The dissipation coefficient becomes 
which always vanishes. 

In striking contrast, non-zero values of the dissipation coefficient can be obtained by using an approximation for 
the X propagator, which we shall refer to as the exponential decay approximation, 

GMsit) « ^-^e-*("'=-^^,)^ ^ > 0^ (42) 

2{ujk - II x) 

where Uk = {k"^ + rrv^Y^'^ and is the x decay width. This approximation gives a non-zero value for the dissipation 
coefficient Tms from eq. (|^ . 

Tms = ^g'cP'lm [ -f^- ^— (43) 

This result was first obtained by Morikawa and Sasaki . 

The difference between to two results is due to the energy dependence of the width and the consequent failure 
of the exponential decay approximation for the zero-temperature propagator. The simplest example is when the x 
propagator has a self-energy contribution Ft from a massless boson loop, which gives 

n = ^ ln(fc2 -uj^- ie) (44) 

where 7 is a constant. The Feynman green function 

Gf{lo) = {-oj^ +ujl + n)-^ (45) 

has branch cuts above and below the real lo axis which push the particle pole off the first Riemann sheet. After some 
contour manipulation, 

/"OO 7 

Gf(<)= / {i-uj^ + Ljl-2irxUJky^ ~ i~uj^ +Lol + 2trxUJky^}e''^' t>0 (46) 

Jk ^7'' 

where F^^, = 7^/(2wfe) and the real part of the vacuum polarisation has been dropped. The result can also be obtained 
directly from the spectral representation of the Wightman function. 
The energy integral can be done analytically, 

Gpit) = ' e-'i'^'^-^^^)^ + ficuk - tTx, t) fiojk + tVx, t) (47) 

Z[UJk — it x) 
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with an error function part 

Aujt ^—iujt 

f{u;,t) = —Eiii{k + Lu)t)-- EMk-^)t) (48) 

We recognise the first term in (|47|) as the earher approximation (|42ll . This term is a good approximation to the full 
expression for early times up to i ~ r^, where 

in 

^i-r-Mog-^, (49) 

X 



when the error functions begin to dominate and the propagator decays as a power law. When 1)47(1 is substituted into 
the formula for the dissipation coefficient H52I) . we find 

r = 4<?V'lm j ^ dt'GF{t')H' = 0, (50) 

due to a cancellation between the exponential and the error function terms. 

A similar argument shows that the equilibrium propagator also has a power law decay at large time. However, 
it may be possible that an exponential decay law is the best representation the non-equilibrium propagators in the 
dissipation term, as has been suggested by Lawrie 0,^1] . There is some support for this suggestion from numerical 
evidence that the real-time spectral function in two dimensions has an exponential decay which extends well beyond 
the relaxation time p3l| . The implication would be that Vms rnay be applicable if the inflaton evolution timescale 
Tcj, lies between the relaxation time and the thermalisation time of the heavy field. This issue awaits further 
investigation. 



D. Non-local dissipation terms 

The zero-temperature propagator has a well-defined relaxation time — T^^, but fails to produce a local dissipation 
term in the inflaton equation due to the non-exponential terms. We shall now investigate the form of the non-local 
dissipation term. We begin by introducing a transform of the inflaton field 

/•OO 

0(c^) = / dtV(t')e'"^*'"*^ (51) 



This substitution gives the linear dissipation term in the equation of motion, 

(Pk 



^^r-'"^ = - Wlm / ^^°°^^^p,(c.l)Px(-2)^(-l+c.,-c.3-^6)-V(-3) (52) 
A more suggestive way of rewriting H52|l is to set 

pnon-loc ^ 2g2jj^ {g'' C^'' K [dt)c^) . (53) 

where the kernel 

K{z) = -lm ——^ ——p^{uJi)p^{uj2)i{uJi+LU2 + iz + ie} (54) 
J [^T^r Jo 27r Ztt 

This would be a local derivative expansion if K{z) was an analytic function. Although K{z) is finite at z = 0, there 
is a branch cut and K{z) is non-analytic in the neighbourhood of the origin. 

We can isolate a leading order logarithmic term in K{z) by taking six derivatives with respect to z and then taking 
a small z approximation. For small z, 

a«A-(.) ^ -120.,(0)^ I ^ |1 ^~ ifi + „,-a ^ ^,,(„,'=-. 

where p^{0) « •^T^/ml. (At zero temperature, vanishes for uj < k, and it is important to take the limit uj k 
and then take the limit k ^ Q). Hence, 

1 

2lmK{z) Ko + K2Z^ + KiZ^ - — ^z-^ln^: (56) 
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where Kq . . . K4 are constants. The important term for dissipation is the logarithm, which appears the field equations 
as 

<^ + 3F<^ + y^-|^^ff2^2 5fln(9O0 = O (57) 

This dissipative term represents the leading order adiabatic behaviour of the full expression (|52|l . 

The non-local operator ln(i9t) is fixed by linearity and causality as discussed in |4J|. It can be expressed in integral 
form, 

In (dt) / = - /* In (e^t - t')) dt'f{t')dt' (58) 

where 7 is Euler's constant. The integral can be evaluated for a broad class of functions, for example 

In (m^^dt) 9{t) sin(m<) = S\{mt) cos(m<) — Gi{mt) sin(r7ii) (59) 

where Si and Ci are sine and cosine integrals. (This can be used to appl y eg . H57fl to the decay of an oscillating field 
(j) = ylsin(mt), where the amplitude is a slowly varying function of time |45|.) 

IV. SELF-ENERGY 

We shall examine the dissipative terms in the equations of motion for the inflaton with a range of interactions which 
are designed to be embedded easily in supersymmetric models. In this section we shall calculate the self energy for 
bosonic and fermionic massive x fields interacting with light bosonic and fermionic fields of mass my. 

The specific interaction terms for the 0{h?) contributions to the self-energy are 



1. For X ^ 2?/ 



2. For X ^ 2y 



CY = -^gH{v\^ +y^\) (60) 



i^Y = --^hixi^yPLlpy + xHyPRi^v) (61) 

3. For x^yy 

Cy = -V2h{y^^PL^y + y^^Pni^y) (62) 
where Pl = ^(1 -I- 75). The fermions used in the supersymmetric models are all Majorana. 

A. Interaction 1: x ^ 2y 

The leading order contribution to the self-energy from the light boson loop is given by figure 01 



P ■ 
X 



a ^ \ b 



y 



FIG. 4; Order contribution to the boson self-energy Ilab. 
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We shall calculate the dissipative part of the self-energy, and denote by tty the function defined in equation ljA5|l . 

ay = i{e^'^ - l)ni2 (63) 

The loop diagram gives, 



-g'h^cb^e^- 1) / ^GMK)GMP " K) (64) 



A spectral representation ljA9p can now be applied, 

ay^lh\e^--l) J |!^^n(fc>(n5(c.-fc"-np(k,fcVk-p,fc°') (65) 

We can use the free bosonic spectral function given in the appendix HA16|1 . The integral reduces to a different form 
for different ranges of energy. Let 



(4m^+p2^i/2. (66) 



In the uj > LUc case, the integral reduces down to 



a,, 



y = 2T:g^h''(b''ief^'^ - 1) / '.^ \' "'"^^"^ Siu; - uou - c^p-fe). (67) 



d^fc 1 + 2n{L0k) 



This integral can be done using a substitution 



y = Wfc + Wp_fc (68) 

X 



ujk-ujp^k- (69) 
The range of integration becomes y > coc and ^x{P) < x < x{P), where 



2\l/2 



In the end, we get 



Repeating the analysis for < lu < lOc gives a = for p < u < Uc and 



/i2 y /I _ e-'3(^(^)+")/2^ 
«y = 7Z5'0'-ln^— ^TTTT^CT^) 0<^<p. (72) 



47r p _ e-/5(^(^')-")/2 

These are quite a complicated expressions, but they can be interpreted simply in certain limits. At zero temperature 
only the first term from the u > Uc result survives. Since a becomes the imaginary part of the vacuum polarisation, 
this term is related to the particle production probability and ujc is the energy threshold for particle production. 

At the other extrerne, for high temperatures a = 0{h^'rn^). This result can be checked by taking the Hard Thermal 
Loop approximation p35l |. which when applied here predicts that the imaginary part of the scalar self-energy should 
have no 0{T^) contribution. This contrasts with the vacuum polarisation obtained from other types of vertex. Gleiser 
and Ramos |33 considered the scalar vaccum polarisation due to quadratic interactions, and their results suggest that 
h'^X^y'^ vertices would give a = 0{K^T^) in the high temperature limit. We shall restrict ourselves to temperatures 
T < h~^m^ for which the 0{h'^T^) terms can be neglected. 

For light fields we can approximate the result by taking the massless limit rriy — 0, 

/,2 7;2 T /I - p-/3('^+P)/2\ 
ay^—g^^^eiuj-p) + —g^q^^-ln ^ . (73) 

This approximation breaks down at the point oj — p, but remains integrable. 
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k 

FIG. 5: Order h'^ fermionic contribution to the boson self-energy. 

B. Interaction 2: x ^ 2y 

The leading order contribution to the self-energy from a fermion loop is given by figure 

We shall calculate the dissipative part and denote this contribution by ay. The fermion loop diagram gives, 

ay = h\l - J ^tr{Su{K)S,^{K - P)). (74) 
A spectral representation ljA9|l can now be applied. In the oj > ujc case, the integral reduces down to 



a 



y = -2nh^ I ^L^I^iP^ + Ami) ^ ^'^ - -P"^)" (^5) 



(27r)3 AuJkUJp^k 

This integral can be done using the same substitution as before, resulting in 

3/2 

Btt" I" ' -P^ / 4:71 I" ' P2 ; g-/3((^-a;(P))/2 _|_ 2 



a, = -^P^l + —l\ -__p2 1 + ^ (76) 



Repeating the analysis for < w < gives 

ay = P^ 1 + n — :^rT^, — ttt^ < uj < p. 77 

y Att p \ P^ j \^e-'3(^(-P)-'^)/2 + ly ^ ^ ' 

The fermion loop, like the boson loop, gives a ~ 0{h?rri^ in the high temperature limit. 
In the massless niy = limit. 

The momenta have been scaled to display the dependence on the temperature in the low temperature limit. Note 
that, at low temperatures, the vacuum polarisation of the x field is dominated by the bosonic loop contribution. 

C. Interaction 3: x ~* yy 

The leading order contribution to the fermion self-energy is given by figure 



FIG. 6: Order h? contribution to the fermion self-energy Eab. 
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We shall calculate the function a, given by equation (|A6|I . 

fi = ~i(e'^" + l)Si2. (79) 

From the diagram we have, 

a = 2h\eP- + ^G,2iK)SMP - K). (80) 

We proceed as before. For a; > t^c, 

^2;^ (1 - n{ujk) + n{<^p-k)) {-l^iPf. - K^) + my)5{uj - uju - ojp^k)- (81) 
The substitution of the new variables x and y can be used again. The gamma-matrix terms become 

-7^(Pm-^m) = ^'^o-^^i, (82) 
where ctq and ai are defined in HA13|) . The integral gives 



with 



and 



/ 4™2y/2y /i_g-/3(^+.(P))x 



(85) 



+^ (Tj (^Li2(e-^("--(^))/2) -Li2(-e-^('^-^(^»/2) _Li2(e-'3(-+-(^'))/2) +Li2(e"'5("+^(^»/2)^ 
where a::(p) was given in eq. H7U|I and 

00 „ 

Li2(^) = E^ (86) 



is the dilogarithm function. 

It is worth noting at this stage that the relaxation time of the x field does not depend on the (Ji terms. This is 
helpful because dissipation coefficient in the small coupling limit is largely determined by the relaxation time. If we 
use eq. (|A15p for the thermal width F-^, with m — — gcj), we get 

2c.pF^ = %V ( 1 + ^1 (1-^) +%V-(l + ^)ln(-^^-^^ Tpss). (87) 

We are interested mostly in very light fields y, for which wc can use the limit my = when 

2.,F, . !^y^^ei.-p) + ^,^0^^1n [^^^^ 

This is identical to the thermal width of the massive boson fields when both the light boson and light fermion loops 
are included. 
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V. MORE DISSIPATION TERMS 



The calculation of the frictional terms in the inflaton equation of motion now needs to be extended to a more general 
class of models. The first generalisation we do is to take into account the coupling of the light fields to the inflaton 
which arises in eq. Ht)U|) . The other generalisation is to consider dissipation terms which arise due to interactions 
between the inflaton and massive fermionic fields. 

The interactions which we consider are all contained in a supersymmetric theory which has three superfields <&, X 
and Y with superpotential, 

W = - Axy2. (89) 

The scalar components of the superfields are 0, x y respectively. This superpotential can easily be modified to 
produce a hybrid inflationary model by adding additional terms or by adding a Z3-term to the potential, neither of 
which affects the dissipation terms. 

The interaction terms responsible for the inflaton decay are 

1. For (f> ^ X 

Ci = -g'c^'xh; (90) 



2. For 



-C/ = -^50V'xV'x- (91) 



The scalar interaction vertices are given by figure |21 



A. Dissipation with an intermediate boson: (f) 2y and (j) ~* 2y 




FIG. 7: Contribution the r] self-energy. 



Consider the dissipation coefficient given by eq. H28|) with the equilibrium self energy diagram shown in figure 
The dissipation coefficient in this case will be denoted by T{4) lyjiy). ft is given to order by 

r(<^^2j//2y) = lm 1 ^^^r/^ri«''''GL,(c.i)G^,,(c.2)2^*a^,<5(c.i+c.2), (92) 

where the V ahc are vertex factors. The vertex factors now include contributions from the light fields, as shown in 
figure ijSJ). The vacuum loop inside the vertex factor is similar to the leading order correction to the x self-energy, 
which allows to write the vertex factor in terms of (Xy^ as 

V^''\lo,uj^,u:^) = -2ig20c2"'' + (t>-\\ + n)aj,(wi)c2'''' + + n)aj,(cj2)c2"'' + 0(/i^) (93) 

where Ca^"" = Sa''Sa%2n)^S{K + Ki + K2), 1 = 2 and 2 = 1. 
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X 

c 



FIG. 8: Contributions to the rj vertex up to order g h 



The vertex factor can be substituted into eq. (|92|l . We define 

a = 2ReGii = -2ReG22 
and use the spectral representation G21 = i{i + n)p-^. The integral reduces down to 



(94) 



(95) 



The first term recovers our result (|34|l from section IIIB, with an extra factor of two because x is now a complex field. 




0.0006 
0.0005 
0.0004 

cF 

S 0.0003 
0.0002 



0.0001 



full expression 
shall h limit 
small T limit 




0.02 0.04 0.06 0.08 0.1 0.12 
T/m 

FIG. 9: The dissipation coefficient T{<j> 2y/2y) plotted as a function of temperature for h = 0.79. The plot on the right 
shows the low temperature region where r{(f> —* 2y/2y) oc T'^ and the small h approximation breaks down. 



Numerical values for the dissipation coefficient are plotted in figure El The small h approximation works reasonably 
well even for /i ~ 1, except in the low temperature regime. 

For small coupling, h, the integrand is dominated by poles in the spectral function. The last two terms in eq. H95|) 
are higher order in h and can be dropped. The first term gives 



r((/. 2y/2y) « g^<l>' f3 



(27r)3c.2 



n{n + 1). 



(96) 
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where is the relaxation time of the x field. For h ^m^ > T, we can use eqs. (|73ll and (|78|) for the boson and 
fermion loop contributions to the self-energy, 

r = '^-^ (97) 

^ ay(ujk,k) + ay{ujk,k) 

When this is combined with a high temperature approximation, the dissipation coefficient becomes 

r((/) -> 2y/2y) w 0.691 1^ T h-^m^ > T :^ m^. (98) 

This is shown in figure El 

In the low temperature limit the integral (|95|l is dominated by small k and lo. We can approximate ~ 2a/'^'* 
and a « The leading order behaviour is determined by the self-energy a^, 

ml J 27r (27r)3 ^ 

We can use the value of ay given in eq. (|75|l . to get 



r(0 -> 2y/2y) « ^^-^ / — -— ^ 34/34a^ n' . (99) 

TO?, 7 27r (ZTT)-^ " 



4 jiS 

r(0 2y/2j/) « 4.0 X 10-2 g2;^4 | j T-^m^. (100) 



TO? 



The value of the coefficient has been found by numerical integration. 

B. Dissipation from an intermediate fermion: (f) ^ yy 

The contribution to th dissipation coefficient from the decay (j) ^ x yy la given by 



r(0 ^yy)=gHm I I dt' tT{S^,{t')Sl{t'))t' (101) 



(2^)^ Jo 

where the fermion two-point functions and conjugates where defined following eq. (|13|l . Using the spectral represen- 
tation of the two-point functions in eq. I|A9|I gives, 

^(^-yy)-2^'jj^l ^^--(^x)- (102) 

In the small coupling limit we can expand the integrand about the poles of the spectral function using (|A15|) . 

r(</. ^ yy) « g^mp f ^ nM{l - n^) . (103) 

The relaxation time = l/F^, where is given by eq. (|88|l . The main difference between this result and the 
contribution from the bosonic decay channel comes from the fermion distribution functions. 

Figure 1101 shows a comparison between the intermediate bosonic and fermionic decays of the inflaton. The total 
dissipation coefficient from all decays at high temperature is 

rw0.97|jT /i-^TOx>T>TOx. (104) 

When T > h^^m^, the two-loop contribution to the x self-energy corresponding to decays such a,s (p ^ vyy becomes 
important. The size of this contribution to the self-energy is known from previous work to be 0{g^h^T'^) |l3| . leading 
to the dissipation coefficient of Hosoya and Sakagami F = 0{h~^m'^T~^) [l^. 

In the low temperature limit the integrals are dominated by small k and w. We have ~ 2a /m<^, and 

r(,^,S,.£l!/|^,^"^8,5W)e^»- (105) 



4 



Using the expression for a in (|84ll and (|85|) . 



F(0 yy) cc g^h^ ^ (106) 



The bosonic decay rate, which is oc T^, dominates at low temperatures. 
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FIG. 10: The bosonic and fermionic contributions to the dissipation coefficient plotted as a function of temperature for h — 0.79. 



C. Dissipation with the exponential decay propagator 

The results have so far assumed that the heavy field has a thermal propagator. We consider now what happens if 
we use the exponential decay approximation for the propagator to calculated the friction coefficient Tms which was 
discussed in section IIIC. Although the heavy field is now not in equilibrium, the light field thcrmaliscs on a different 
timescale which we suppose is sufficiently short to keep the light field in equilibrium. 

For the bosonic decay we obtain the contribution Tms{4> ~^ '^v) to the dissipation coefficient from the decay 
X 1/2/ by using eq (O, 

rA/s(<^ - 2y) = 16ffV / (0^^' (107) 
The self-energy Uy is determined by the results of section IV. The dissipation coefficient can be obtained analytically, 



where 

2 

^ Jo (1 + a2a;2)4 - 1 
Similarly, the contribution rj\/5 {4> — > 2y) to the dissipation coefficient from (j) ^ X ^ VV becomes 



/(") = - I „ , ..2^2., .. , ■ (109) 



lOTT^ V mv / \ V mv / V 



where 



2 

TT Jo (1 + a'^x'^Y + 1' 



/»-z/ n , .2..,4.. , (Ill) 



The main application of these results would be at low temperature and the temperature corrections are relatively 
small. 
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For the fermionic decay 4> ^ X ^ UV^ can make use of the symmetry between the decay width of the fermion x 
given by eq. I|88|) and the decay width of the boson x when = to^ . We obtain the contribution to the dissipation 
coefficient, 



rMs(0-yy) = ^™^ 1+ — / — (112) 




In the zero-temperature hmit, the total dissipation coefficient is then 

Ta/s = 2 + V — 5- g(j). (113 

This is similar to the results obtained for non-supersymmctric models by Morikawa and Sasaki 8] and by Berera 
and Ramos js^l , although the coefficient calculated here is significantly larger due to the large number of decay 
channels available to inflaton in the supersymmetric theory. In addition, there are small numerical differences between 
different results due to slightly different prefactors in the exponential decay approximation. 



VI. ENTROPY PRODUCTION DURING INFLATION 



The main application for our results is to the inflationary universe. In this section we shall make some general 
remarks and check the consistency of some of the assumptions that we have made. During the slow evolution which 
is characteristic of inflation, we replace the inflaton equation by 

3i7(l + r)0 + l/' = O (114) 

We have parameterised the effectiveness of dissipation by taking the ratio of the dissipation coefficient to the expansion 
rate, 

r=^. (115) 

The energy density during inflation is dominated by the inflaton potential energy V ^ and the expansion rate is given 
by 

3i72 ^ (116) 



By combining these equations, we can compare cj) to H, 



A2 g ^2 



p 



47r(l + r) 

where e measures the slope of the potential, 



(117) 



ml fV'V 

The second time derivative of the scale factor is positive for e < 1 + r. 

The temperature of the radiation field can be determined from energy momentum conservation. The energy lost 
by the inflaton is converted into heat and increases the entropy density of the radiation s 34] , 

Ts + 2,HTs^T4? (119) 

The radiation is a relativistic gas, with 

s = ^aT^ (120) 

where a is the effective radiation constant, taking into account the number of degrees of freedom. We see immediately 
that for r oc T^, the temperature does not fall to zero, but reaches a lower limit T^, 

TrS = r^^. (121) 
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If r falls off faster than T^, then the cosmological redshift would drive the temperature downwards and the radiation 
into the vacuum state. 

We shall consider this low temperature regime now in more detail, using the thermal expression for the dissipation 
coefficient, 



^3 



ig'h'^ — . (122) 
ml 



The temperature falls to the value given by eq. (|121|l . When we eliminate using (|117|l . 



2 



H Wna (l + r)^TO?, 



(123) 



Further information can be obtained from the amplitude of fluctuations in the cosmic microwave background. These 
have their origin as thermal fluctuations during inflation ii T > H. Theory gives an amplitude A for a given mode, 
where 

evaluated at the time that the mode crosses the horizon, and c « 0.045. The measured amplitude A « 5.4 x 10^^ on 
a 500 MPc scale relevant to the large scale structure. We can use A to eliminate the expansion rate and obtain 

TTip \256tt ac^ ) 1 + r ^ 

The range of temperatures for which the dissipation coefficient calculation can be applied and still be consistent 
with the low temperature approximation is H < T < m^. This translates to limits on the values of which are 
consisten when 

gh^ > (^^^ ^-i/2a2/3 « 5.1 X 10^2 (126) 
It is quite possible for supersymmetric models of inflation to have coupling constants in this range. 



VII. CONCLUSIONS 



We have taken a detailed look at the linear dissipation term in the equation of motion for an inflaton field embedded 
in a supersymmetric model. The dissipation was caused by the decay of the inflaton into light particles during inflation 
with a heavy field x acting as intermediary. 

• For temperature T > m^ , our results are consistent with a wide body of previous work. Some differences in the 
temperature dependence can be attributed to the special features of the interactions which play an important 
role for temperatures close to the heavy particle mass. 

• For H < T < m^, where H is the expansion rate, we find significant differences between the dissipation 
coefficient in the thermal background and the dissipation coefficient obtained assuming exponential real-time 
decay of the x propagator [MEMI^- In particular, the dissipation coefficient in the thermal background does 
not approach a constant as T ^ and H ~^ 0, but instead behaves as a power law cx T^. 

Thermal propagators decay exponentially for only a limited time before switching over to a power law decay. The 
switch over occurs at a logarithmic factor times the exponential decay time. A purely exponential decay might 
possibly represent an approximation to the non-equilibrium behaviour of the system if the dynamical timescale of the 
inflaton where shorter than the thermalisation timescale of the heavy field. We have therefore given results for both 
the thermal and exponential propagators. 

We have argued that the low temperature dissipation coefficient obtained using the thermal propagator can still 
have an important effect on inflationary models. For a signiflcant range of parameters, the dissipation heats the 
radiation to the H < T < m^ range. Thermal fluctuations are then more important than vacuum fluctuations and 
become the source of primordial density perturbations P, |H IE ^i 1^ QE 03 ■ 
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Dissipation is far stronger for the exponential propagator and warm inflation is realised for a far larger range of 
parameters. This case has been studied in previous work which has divided up the parameter ranges into those 
leading to cold and warm inflation |23, l3ll Is^ . There remains much to understand about the non-equilibrium physics 
of reheating during the inflationary era, but the results which we have obtained so far indicate that dissipation is 
important during inflation for supersymmetric models with moderate values of the coupling constants. 
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APPENDIX A: PROPERTIES OF THERMAL PROPAGATORS 

This appendix lists some important properties of the thermal equilibrium propagators. We begin with some of the 
most basic rules, 



Kubo-Martin-Schwinger periodicity relations. 



n2i = e'^'^Hu (Al) 
S21 - ~e^'^^i2 (A2) 



• The Kobes-Semenoff rules. 



im(nii) - -^(ni2 + n2i) (A3) 

MSn) = -^(Si2 + S2i) (A4) 

The Kobes-Semenoff rules follow from the maximum time rule. 

The periodicity relations imply that the off-diagonal components of the self-energies can be expressed in terms of 
functions a{P) and a{P), 

1121 ~ + ni2 = ~ina (A5) 

S21 = — «(1 — n)a S12 — ina (-A-6) 

where 

n(uj) — —-. -, ~ —5 (A7) 

are occupation numbers for bosons and fermions. The off-diagonal components of the full two-point functions can 
also be expressed in terms of a real function, the spectral density p, 

G2i^i{l + n)p Gi2=inp (A8) 
S'21 = j(l - n)p S12 = -inp. (A9) 

If we discard the real part of the vacuum polarisation, then the Kobes-Semenoff rules imply that the spectral density 
p and the function a are related by 

The spectral density typically has a rich variety of poles and branch cuts. For a pole at w = ± ir^^, we can 
regard uip as the energy of the particle state and r as the relaxation time. In perturbation theory, ujp — {p^ + m^)^/^ 
and T — 2ijjpa{p,ujp)~^ . 

For fermions, the spectral density 

p = i{-^^Pf^ -I- TO -t- iay^ - ii'j'^Pfj^ + m - (All) 
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We can use rotational invariance to express a in terms of three scalar components, 

a = OLQao + aiCFi + amtn (-^-12) 

where the gamma-matrices have been adapted to a momentum frame, 

o-Q = 7°w - 7 • pp (A13) 
(Ti = 7-pa; — (-^-14) 

These combinations satisfy —a^ = erf = and {ccfTi} = 0. The poles of the spectral function are again given 
approximately by ±Up ± ir"^, but now 

WpT"^ = (ao + Q!m)m^ (A15) 

to leading order in a. The relaxation time does not depend on ai. 

Spectral densities for a free theory can easily be obtained by taking the limit a ^ 0, 

p = 27rsgn(w)(5(p2 +m2) (A16) 
p = 27rsgn(a;)(-7'^P^ +TO)(5(P2 + m2). (A17) 

The limit depends on the sign of co because a{p, —w) = —a{p,u). 
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